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Structure properties of 226 Th and 256 > 258 > 260 F m fl ss j on fragments: 
mean field analysis with the Gogny force 
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The constrained Hartree-Fock-Bogoliubov method is used with the Gogny interaction D1S to 
calculate potential energy surfaces of fissioning nuclei 226 Th and 256 > 258 . 260 F m U p to very large 
deformations. The constraints employed are the mass quadrupole and octupole moments. In this 
subspace of collective coordinates, many scission configurations are identified ranging from symme- 
tric to highly asymmetric fragmentations. Corresponding fragment properties at scission are derived 
yielding fragment deformations, deformation energies, energy partitioning, neutron binding energies 
at scission, neutron multiplicities, charge polarization and total fragment kinetic energies. 

PACS numbers: 21.60.Jz, 24.75.+i, 27.90.+b 



I. INTRODUCTION 

Our knowledge of the fission process has made huge 
progress in recent years with the measurement of mass 
and charge distributions of fission fragments for 70 fis- 
sioning systems [![, performed at the secondary beam 
facility at GSI. The measured fragment yield distribu- 
tions have revealed new kinds of systematics on shell 
structure in nuclear fission, such as transitions from 
single- and double-humped mass distributions to a triple- 
humped structure in the vicinity of 227 Th. From a theo- 
retical point of view, microscopic self-consistent meth- 
ods appear to be well suited to study structure ef- 
fects in fissioning systems, where the sole input is the 
nucleon-nucleon force. Many studies based on mean 
field approaches using Gogny or Skyrme forces have re- 
cently been devoted to the different fission modes, as 
for example in 256-258 Fm isotopes [| g [| [| g [7^ 

where bimodal fission has been experimentally identi- 
fied 0, Jl M El [H [H, [H, El EE E3 and an- 
alyzed [la . Il9| . Furthermore, two-dimensional time- 
dependent calculations have also been performed for the 
238 U isotope in the elongation-asymmetry plane, where 
it appears that fragment mass and total kinetic energy 
distributions are well reproduced. These calculations 
have employed the Time-Dependent Generator Coordi- 
nate Method treated at the Gaussian Overlap Approxi- 
mation and used Hartrcc-Fock-Bogoliubov states [2(| • 

The present work, based on the constrained Hartree- 
Fock-Bogoliubov (HFB) method and the D1S force, is fo- 
cused on the calculation of structure properties of nascent 
fission fragments of light and heavy actinides, namely 
226 Th and 2 56,258,260p m _ Fragment deformations, de- 
formation energies, energy partitioning, neutron binding 
energies, neutron multiplicities, charge polarization, and 
total fragment kinetic energies are calculated for a wide 
range of fragmentations. This large scale study has been 
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made possible thanks to the new generation of fast com- 
puters made available to our laboratory. By the mean- 
time, it is hoped that the calculated structure informa- 
tion here collected for a wide variety of fission fragments 
will serve as guideline for updating inputs (excitation en- 
ergy, energy partitioning, neutron binding energy, etc. . . ) 
to phcnomenological evaporation models aimed at calcu- 
lating prompt neutron emission from, and 7-ray decay of 
fission fragments [H [22|, [23[ . 

The paper is organized as follows. In Sec. II is outlined 
the constrained HFB method in which, like in Ref. [20l ]. 
quadrupole and octupole mass operators are adopted for 
external fields. This section also presents the mean field 
methods used to describe: i) the scission mechanism as 
well as ii) nascent fission fragments in low energy fis- 
sion. In Sec. Ill results are discussed, among which po- 
tential energy landscapes, scission configurations, and fis- 
sion fragment properties. Fission fragment yields arc not 
considered in this work as they require a dynamical treat- 
ment [2^] . Comparisons are made between present pre- 
dictions and experimental data for total fragment kinetic 
energy ( 226 Th, 256 Fm) and prompt neutron multiplicity 
( 256 Fm) of fission fragments. 



II. SELF-CONSISTENT APPROACH TO 
SCISSION 

A. Constrained Hartree-Fock-Bogoliubov method 

The deformed states of the nuclei under study have 
been determined using the constrained Hartree-Fock- 
Bogoliubov (HFB) [24| theory based on the minimization 
principle of the energy functional, namely 

S < <5>({qio})\H-\NN-\ z Z-J2^Qio\H{qio}) >= 0, 

where H is the nuclear microscopic Hamiltonian, Qiq a 
multipolc operator, and \n, \z-, an d A; the Lagrange 
parameters associated to constraints on nucleon numbers 
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N, Z and average deformations qio, respectively, 

< mq la })\N\<i>({q la }) > =N, 

< ${{q l0 })\Z\*({q l0 }) > = Z, (2) 

< H{qm})\Qio\$>({qio}) > = oto. 

and where Q;o is defined as 

r~s — A 

4o = (i + My2iTiE r ^^'^)- ( 3 ) 

z— 1 

In the present study, the Hamiltonian _ff is built using the 
finite ran ge a nd density-dependent nucleon- nucleoli D1S 
force [25l |26||. One-body and two-body corrections for 
center of mass motion are taken into account in H . Many 
calculations have shown that the energy functional de- 
rived from this Hamiltonian provides a very satisfactory 
reproduction of nuclear properties over the whole mass 
table [I?) and especially in the actinide region (28[. In 
Eq. |(3J) the set of constraints {Qm} includes the isoscalar 
axial dipole, quadrupole, and octupole mass moments 
Q10, Q20 and Q30, respectively. The dipole moment has 
been constrained to zero so that the mean position of the 
nucleus center of mass is located at the origin of the coor- 
dinate system. The HFB energy of the deformed system 
is defined as 

-Ehfb(920? 930) = < *(<?20, 93o)|^|*(g20, <73o) > ■ (4) 

In the present study, the Bogoliubov space has been re- 
stricted by enforcing axial symmetry along the z— axis 
and the self-consistent Tn 2 symmetry, where T is the 
time-reversal operator and II2 the reflection with respect 
to the xOz plane. The system of Eqs. |T|) and {2} has 
been solved numerically by iterations for each set of de- 
formations by expanding the single particle states onto 
axially symmetric harmonic oscillator (HO) bases. For 
small elongations (520 < 200 b) a one-center HO basis 
with TV = 14 major shells has been used while for large 
elongations ((720 > 200 b) a two-centers HO basis with 
twice N = 11 major shells has been preferred [2^|. The 
parameters of the one- and two-centers HO bases have 
been optimized for each set of deformations, and we have 
checked that the basis sizes are large enough. The poten- 
tial energies discussed below are the HFB energies defined 
in Eq. Q. 

B. Scission mechanism in the (920 > 930) plane 

At large quadrupole moment, it becomes energetically 
more favorable for a fissioning system to split into two 
separated fragments, rather than to take on a very elon- 
gated shape with a neck. In deformation space, this tran- 
sition corresponds to an evolution from the so-called fis- 
sion valley to the so-called fusion valley [3(| . If a point A 
in the fission valley leads to a point B in the fusion valley 
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FIG. 1: Comparison between symmetric fission of the Th 
and 256 Fm nuclei through total energy, hexadecapole moment 
and minimal density along z— axis in the neck pn- 920 an d 
Z?' ps ^ represent elongation and HFB energy of the first post- 
scission (ps) point for each fissioning system, respectively. 



through a small increment in either one of the deforma- 
tion parameters, then point A is here defined as a scission 
point and point B as a post-scission point. Unfortunately, 
there is no universal way to distinguish a point in the fis- 
sion valley from a point in the fusion valley, and several 
criteria have been used in previous studies to achieve this 
classification. For instance, Bonneau et al. [3l|, H2] con- 
sider that scission occurs when the nuclear interaction 
between fragments is less than 1 % of the Coulomb repul- 
sion energy, whereas in refs. [10, H(| it was noted that the 
scission mecanism in 238 U and 240 Pu is associated to the 
following three properties: i) the neck between the frag- 
ments suddenly vanishes, ii) the hexadecapole moment of 
the system decreases, and hi) there is a drop in the poten- 
tial energy of the fissioning system. Whereas these three 
criteria appear to be equivalent in the U-Pu region, this 
is no longer the case for some of the nuclides studied here. 
As an example, in Fig. [T]wc show the evolution of HFB 
energy, mass hexadecapole moment and minimal density 
in the neck as functions of the quadrupole moment for 
the symmetric fission of 226 Th and 256 Fm. In the top 
panel, the evolution of the HFB energy shows that scis- 
sion can either correspond to a sudden loss ( 226 Th curve) 
or to a smooth decrease ( 256 Fm curve) of the binding en- 
ergy. The same difference in behaviors can be observed 
for the evolution of the mean values of the hexadecapole 
moment (Q40). These examples clearly illustrate that 
scission points can be determined neither by an encrgy- 
nor an hexadecapole moment-based criterion for the scis- 
sioning system in 256 Fm. The lower panel of Fig. Q] shows 
that in the symmetric fission of both 226 Th and 256 Fm, 
the density in the neck displays two different values be- 
fore and after scission, with an abrupt drop at scission. 
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In this situation, we define a post-scission configuration 
as one for which in the matter density along the sym- 
metry axis there is a local minimum that is lower than 
p = 0.06 fm -3 . Using this criterion, we define for each 
nucleus a set of scission points in the (<?20; Q30) plane, 
that is called the scission line. 

Depending on the nucleus and fragmentation, the scis- 
sion transition is cither smooth (e.g. symmetric frag- 
mentation of Fcrmium isotopes) or abrupt, in the present 
subspace of collective coordinates. In the first case, the 
energy of the fissioning nucleus evolves smoothly without 
discontinuity from outer saddle to scission, becoming the 
Coulomb repulsion between nascent fragments at large 
elongation. In the second case, there is an abrupt de- 
crease of energy and hexadecapolc moment at scission. 

Figures and show the evolution of the nuclear 
density at large elongation for the symmetric fragmen- 
tation of 226 Th and 256 Fm, respectively. While 226 Th 
displays a very elongated shape at the scission point (up- 
per panel) and two prolate fragments at the post-scission 
point (lower panel), 256 Fm symmetric fission leads to 
two nearly spherical fragments separating smoothly. In 
the literature, these quite different ways of fissioning 
are called Elongated Fission (EF) and Compact Fis- 
sion (CF), respectively [H, In the 256 Fm sym- 
metric fission case, CF is currently explained by the 
proximity of double-magicity of the fragments (Z = 50, 

n = 82) dun. 
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FIG. 3: Same as Fig. [2]for the symmetric scission of 2j6 Fm. 



C. From scission to fragments 

The main purpose of identifying the scission configu- 
rations of the nuclear shape using the method described 
in the preceding section is to obtain information on frag- 
ment properties and distributions. The underlying as- 
sumption is that, once a scission configuration is reached, 
splitting of the nucleus will occur irremediably yielding 
two separated fragments moving away from each other 
under the action of their mutual Coulomb repulsion. Ob- 
servable fragment properties such as kinetic energy or 
excitation energy can then be inferred from the charac- 
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tcristics of the nascent fragments - as distance between 
centers of mass, deformations, ... - at scission. It is im- 
portant to stress that the fragment properties derived in 
such an analysis will not necessarily all correspond to 
those observed in experiments, since some of the config- 
urations found at scission may not occur with significant 
probability in the fission process. 

For each scission point, a sharp cut is made at the neck 
position zn on the z— axis, which serves to define the light 
(L) and heavy (H) fragments. Some fragment properties, 
namely quadrupole and octupolc deformations, masses 
and charges, distances between centers of charge and 
mass, are next calculated as mean values 

/zn roc 
dz / r.dr6p(z,r), (5) 
-co JO 



(d) H =2TT dz r.dr6p{z,r), (6) 

Jz N JO 

where p is the nuclear density and O a one-body operator. 

We have checked that all first multipole moments from 
Q20 to Qqo of the fissioning system are continuous along 
scission lines, which ensures that the scission configura- 
tions analyzed form a continuous set from which fragment 
properties can be consistently derived. At this stage a 
few remarks are in order, namely: i) the adopted sharp 
cut assumption inevitably leads to non-integer values for 
calculated fragment charges and masses; ii) as our model 
is restricted to two collective coordinates, only one scis- 
sion configuration is predicted for any fixed (92s, <?3s) 
value. As a consequence, the set of fragment pairs here 
deduced is only a fraction of all possible pairs which 
would be formed if the constrained HFB calculations 
were extended to include other collective coordinates. 
Nascent fragments associated with different scission con- 
figurations may be found having nearly the same proton 
and neutron numbers. As charge and mass fission frag- 
ment yields are outside the scope of the present static 
model, such fragmentations are considered having the 
same weight in figures shown below where they will dis- 
play multiple- values, for example when plotted as a func- 
tion of fragment mass. 

III. RESULTS 

A. Potential energy landscapes 

The potential energies have been calculated on a (920 7 
tfeo) mesh from (q 20 = 0, q 30 = 0) to (q 20 = <?2s, 930 = 
93 S ), where (q2 a , l3s) belong to the scission lines. With 
the chosen mesh dimensions Aq 2 Q = 10 b and Ag3o = 
4 b 3 / 2 , each of the potential energy landscapes shown in 
Fig. Q] are generated with approximately 600 calculated 
values. For convenience, the range of potential energies 
shown is limited to 20 MeV for 226 Th (see Fig. g^a)) and 



to 50 MeV for the three Fm isotopes (Figs. Sfb), BJc), 
SJd)). Isolines are separated by 1 MeV. 

The topological properties displayed by the four land- 
scapes are quite contrasted. We first notice that the low- 
est potential minima of 226 Th and 256 " 260 Fm are all soft 
against quadrupole and octupole deformations, which 
should favor coupled quadrupole and octupole vibrations 
at low excitation energies. These are the common fea- 
tures expected for these nuclides at normal deformations. 
As axial deformation increases beyond the inner barrier, 
a well defined superdeformed (SD) potential minimum 
is taking place only for 226 Th. The SD potential mini- 
mum is vanishing for 256 - 260 p m as discussed previously 
for actinides with neutron number N > 156 [281 ] . 

Beyond the SD potential minimum, a valley a few MeV 
deep is showing up in 226 Th for asymmetric deformation 
all the way to a scission point with large left/right asym- 
metric fragmentation. An isomeric minimum appears for 
1720 = 140 b, 930 = 20 b 3/2 . At elongation q 2Q > 150 b, a 
symmetric valley is also observed until the scission point 
q 2s = 500 b is reached. As scission energies are similar 
in both valleys, symmetric and asymmetric fission modes 
are expected to compete in this nucleus. For 256 - 260 ;F m , 
the potential landscapes display similar and smooth pat- 
terns beyond the first axial barrier. In contrast to 226 Th, 
we observe that: (i) a shallow asymmetric valley is iden- 
tified for <7 30 > 30 b 3 / 2 , (ii) the fall-off of the poten- 
tial landscapes versus elongation for the latter nuclides 
is smooth for asymmetries 530 > 50 b 3 / 2 , (iii) the scission 
lines display approximately smooth and linear trajecto- 
ries over the (q 2 o, <73o) plane, and (iv) a symmetric valley 
is gradually developing beyond q 2 Q = 100 b as N grows 
from 156 to 160. This last feature is not inconsistent 
with the observation of a transition from asymmetric to 
symmetric mass division in fission, in going from 256 Fm 
to 258 Fm [HIM HI Whether or not this transition 
can be further analyzed with the present static mean field 
approach will be discussed below. 



B. Scission lines over the (520, ?3o) plane 

The mesh sizes A.q 2 $ and Ag3o so far adopted are well 
suited for performing a survey of potential energy land- 
scape properties. In the vicinity of scission points the 
step sizes have been dramatically reduced to Aq 2 Q = 2 b 
and A</3o = 1 b 3 / 2 , in order to define scission points with 
high precision. For each nucleus, approximately two to 
three hundred scission points are used to define a scis- 
sion line. This is illustrated for 226 Th in the upper panel 
of Fig. where each point in the (q 2 o, 930) plane corre- 
sponds to a single HFB calculation. Only configurations 
before scission are shown. The curve in red color is the 
scission line, which is made of all exit points (q 2s , <?3 S ). To 
ease forthcoming discussions, a few scission points have 
been labeled with letters a, b, c,. . . , j. 

Most of the constrained HFB calculations at given (q 2 Q , 
930) values are performed using as a starting point the 
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generalized density matrix 1Z [3^| obtained at {q 2 n — Aq 2 Q, 
qso). However, in a few cases, it has been necessary to 
start from the density matrix calculated at either (520 j 
q 30 - Aq 3a ) or (q 20 + Aq 20 , <73o), hi order to reach all pos- 
sible fragmentations. For example, the segments defined 
between the labels b and c and between the labels d and 
e were determined increasing asymmetry and decreasing 
elongation, respectively. 

The scission lines determined for 226 Th and for 256 Fm, 
258p m anc j 260p m are sn0 wn in the upper and bottom 
panels of Fig. respectively. The lines for the Fm iso- 
topes display similar features. The symmetric scission 
configurations are found at q 2s = 270 b. Beyond this 
point, q 2s and q 3s increase gradually until q 2s reaches a 
maximum for q 2s ~ 500 b where asymmetry takes on val- 
ues in the range q 3s = 80 — 100 b 3 / 2 . For higher asymme- 
tries, the scission lines display wiggling patterns and are 
quite similar. The trajectory followed by the 226 Th scis- 



sion line over the (q 2 o, (730) plane is quite different. First, 
the nucleus stretches and gets an elongation nearly twice 
as large as the one for Fm nuclides before symmetric scis- 
sion takes place. Next, elongation decreases as asymme- 
try increases until q 2s reaches a minimum for q 2s = 250 b, 
(label c in Fig. 0. Except for the point on the scission 
trajectory marked with the label f, q 2s and q 3s increase 
smoothly until the point labeled i is reached. Beyond 
this point located at (q 2s = 444 b, q 3s = 142 b 3 / 2 ), both 
q 2s and q 3s decrease until the scission point labeled j is 
reached. Although the scission line is defined beyond the 
point labeled j, this segment lies in a (520, ?3o) region 
where potential energy is sharply raising. Therefore, the 
corresponding scission configurations will not be reached 
in low-energy fission and they will not be considered in 
the rest of this work. For the same reason, the Fm scis- 
sion points beyond (q 2s = 410 b, q 3s = 111 b 3 / 2 ) will also 
be discarded. 
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FIG. 5: (Color online) Upper panel: the scission line for 226 Th 
is shown over the (920, Q3o) plane as a continuous curve in 
red color along which are marked symbols a, b, c,. . . j. The 
black dots, representing single HFB calculations, are shown to 
illustrate the densening of the mesh used close to the scission 
line. Lower panel: scission lines for 256-260 Fm. 



C. Energy along scission line 

The potential energies -Ehfb along scission lines are 
shown for 226 Th and 256 - 260 F m as functions of the frag- 
ment mass Afrag in Figs. [5] and respectively. These 
energies take on identical values on both sides of the 
symmetric fragmentation where Af rag = A/2. On these 
figures, each solid dot is for a single HFB calculation. 

In 226 Th one principal and two secondary minima are 
observed which arc likely to represent the most probable 
fragmentations in low energy fission. Hence, both sym- 
metric (Af ra g ~ 113) and two asymmetric modes (Af lag ~ 
132 and A[ lag ~ 145) are expected for this nucleus. Frag- 
ment charges in the symmetric mode (Zf rag ~ 45) and 
the two asymmetric modes (Zf rag — 52 and Zf r&g ~ 57) 
appear in good agreement with those found in the triple- 
humped mass/charge distribution measured a few years 
ago and analyzed in terms of the superlong (Zf rag = 45), 
standard I (Zf rae = 54) and standard II (Zf ras = 56) 
fission channels 14011. 
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FIG. 6: Th. Potential energy along the scission line as a 
function of fragment mass. The symbols a, b, c,. . . j have the 
same meaning as in Fig. [5] See text for more details. 



Figs. O and [11] show that there is a correlation between 
the structures in the potential energy along the scission 
line and the behavior of the scission line in the (920 , <?3o) 
plane. In order to better visualize this correspondence, 
the potential energy of characteristic scission configura- 
tions labeled as a, b, c,. . . in Fig. [5] is displayed in Fig. [6] 
One observes that, as Af Tag increases, i) the scission line 
shifts from symmetric to asymmetric mass division fol- 
lowing an irregular trajectory over the ((720 , Q30) plane, 
and ii) to each labeled scission configuration is associated 
a break in the -Ehfb energy values. It thus seems that the 
competition between symmetric and asymmetric fission 
of 226 Th is tied with the static structure properties of the 
fissioning system along the scission line. The asymmet- 
ric scission configurations calculated for Af Tag — 132 and 
Afrag — 145 coincide with the points marked with the 
symbols f and i, respectively, in Fig. [5J 

The absolute minima in -Ehfb for 256 ~ 260 Fm along 
the scission line take place for asymmetric fragmentation 
with Af lag ~ 145, property which correlates rather well 
with the location of a peak in the fragment mass distribu- 
tions identified for Af rag — 142 in the 256 Fm mass-yield 
measurements l_l( . Symmetric fission is not energetically 
favored as -Ehfb displays a maximum for Af rag = A/ 2, in 
contrast to the above results for 226 Th. However we ob- 



serve that the difference in energy between the maximum 
and minimum values taken by -Ehfb for -Afrag = A/2 
and A fr ag ~ 145 decreases from 22 McV to 16 MeV 
as the mass of Fm isotopes increases from A = 156 to 
A = 160. Although this feature would favor a transition 
from asymmetric to symmetric fission, making a more 
definite conclusion on this transition requires a full dy- 
namical calculation in which both potential energy and 
tensor of inertia from ground state deformation to scis- 
sion configurations play a role. 
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FIG. 8: (Color online) Axial mass quadrupole moments (Q20) 
of the nascent fission fragments for 226 Th and 256 " 260 Fm. 



D. Fragment deformations 
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FIG. 9: 112 Ru and 150 Ce potential energy curves from con- 
strained HFB calculations restricted to axially-symmetric and 
left-right-symmetric shapes as a function of axial quadrupole 
deformation [4l| . The potential energy of 150 Ce has been ar- 
bitrarily increased by 295 MeV to ease comparison between 
curves. 

formation corresponds to a shallow secondary minimum 
of the potential energy curve of the Ruthenium isotopes 
as a function of quadrupole deformation, as illustrated in 
Fig. [9] for 112 Ru. Very heavy fragments around A ~ 170 
are also predicted to be well-deformed. This shell effect 
is associated to the deformed magic numbers Z = 58 and 
N = 92 at (Q20} — 15 b. The potential energy curve for 
1 5gCe is also plotted in Fig. [9] as a function of the axial 
quadrupole moment, and it appears that {Q20) — 15 b 
corresponds to the ground-state deformation. 

Axial mass octupole moments of fission fragments are 
plotted in Fig. [10] as functions of the fragment mass. The 
octupole moments display almost the same behavior ver- 
sus -^frag the one for the quadrupole moments: minima 
are observed for A ~ 86 and A ~ 130 and maxima for 
A ~ 112 and A ~ 170. 



The axial mass quadrupole moment of the nascent fis- 
sion fragments along scission lines is plotted on Fig. [5] 
for the four studied fissioning systems. The most strik- 
ing feature is that the fragment deformations do not 
significantly depend on the fissioning system. The four 
curves are almost superimposed and have the expected 
saw-tooth structure: minima are found for A ~ 86 and 
A ~ 130, and maxima for A ~ 112 and A ~ 170. In- 
deed, on the one hand, strong spherical shell effects for 
N = 80 and Z = 50 stabilize spherical fragments of 
Tin isotopes at scission. In the case of 226 Th a mini- 
mum with (Q20) — 5 b is also observed, corresponding 
to Krypton isotopes with A ~ 86. This effect is driven 
by the neutron magic number N = 50. On the other 
hand, well-deformed Ruthenium isotopes (Z = 44 and 
N ~ 68) are here predicted with (Q20) — 22 b. This de- 



E. Fragment deformation energy 

Energy partitioning in fission is a key input of models 
aiming at describing sequential neutron and 7-ray emis- 
sion from fission fragments [111, [13, H3| . In the present 
study, the assumption will be made that the excitation 
energy stored into fission fragment arises only from their 
quadrupole and octupole deformations at the moment 
of scission. With this assumption possible intrinsic or 
thermal excitations prior to scission are neglected. The 
estimates given below must therefore be considered as 
lower bound of fragment excitation energies. 

The fragment deformation energy is defined as [IH 

Edci = Eg — E gs , (7) 
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FIG. 10: (Color online) Axial mass octupole moments (Q30) 
of nascent fission fragments for 226 Th and 256_260 Fm. 



where Eg is the energy of the nascent fragment, and E gs 
the one of the fragment ground state. In this work, E gs 
has been deduced for all fragments from usual HFB cal- 
culations, whereas Eg is the HFB energy predicted in a 
constrained HFB calculation where the axial quadrupolc 
and octupole moments are those obtained at scission con- 
figurations (see Figs. [81 and [T0|) . Let us mention that in 
these two sets of calculations, the neutron and proton 
numbers of each fragment have been taken to be integer 
values closest to the N and Z mean values calculated for 
the nascent fragments. Such an approximation leads to 
an uncertainty in HFB energies which amounts to be less 
than 1 MeV. 

The FF deformation energies (i?def) derived in this way 
for the four nuclei studied here are shown as functions of 
Af rag and Zf reig in Figs. fTlTa) and [TTTb). respectively. 
Strong variations are observed, with maxima reaching 
^def ~ 15 — 20 MeV near Af rag ^120 and minima close 
to zero near Af Yag ~ 145 and A[ rag ~ 130. This latter 
minimum corresponds to symmetric division in Fm nu- 
clides. Its origin is of course the occurrence of strong 
shell effects in nuclei close to 132 Sn. 

When plotted as function of Zf Tagl regions with £?de£ ~ 
correspond to Zf Ta , g ~ 50 and Zf Ta , g ~ 56, that is 
to 128 - 130 S n ail d 144 Ba, respectively. Furthermore, the 
maximum identified previously in the E^d values at 
^4frag ~ 120 gets split over two Zf rag components, namely 
Zfr ag ~ 48 and Zf rag ~ 52, that is for near symmetric and 
highly asymmetric charge divisions in the Fm and Th nu- 
clides, respectively. 

Finally, Fig. HUc) displays the difference (E^ - En) 
between the deformation energies of light and heavy 
fragments. This difference takes on values ranging 
from 23 MeV to -15 MeV. Extrema are located at far- 
asymmetric mass divisions. More than 70% of the light 
fragments display E^ > En values. 
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FIG. 11: (Color online) Nascent fragment deformation ener- 
gies for 226 Th and 256 ~ 260 Fm as functions of: (a) fragment 
mass, and (b) fragment charge. The differences between light 
(L) and heavy (H) fragment energies as functions of light frag- 
ment mass are shown in (c). 
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FIG. 12: (Color online) One-neutron binding energies of 
nascent fission fragments as functions of fragment mass for 



'Th and 



'Fm. 
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FIG. 13: (Color online) Differences between one-neutron 
binding energies of fragments as calculated for scission con- 
figurations and for ground states, plotted as functions of frag- 
ment mass. 



F. Prompt fission neutrons 

In the present section, we aim at calculating the mul- 
tiplicity ^f rag of prompt neutrons emitted by each fission 
fragment. For this purpose, we assume that the defor- 
mation energy of any fragment is converted into internal 
excitation energy through collective vibrations and that 
the fragment will de-excite only through prompt neutron 
emission. As an estimate, the neutron emission multi- 
plicity of one fragment is taken as (U [3] 



where B* is the one-neutron binding energy in nascent 
fragment, and (Ek) the mean energy of the emitted neu- 
tron. The latter is assumed to be 2 MeV in 226 Th 
and 1.5 MeV in 256-260p m ggj 

1. One-neutron binding energy 

In the present work the one-neutron FF binding energy 
B* is taken equal to the neutron chemical potential ob- 
tained in the HFB calculations performed on the scission 
line. The B* values arc plotted in Fig. [T2] as a function 
of fragment mass. It is seen that they globally decrease 
from approximately 7 MeV to 3 - 4 MeV with increasing 
mass. The lowest B* values are obtained for A ~ 136 
(Z ~ 50, N ~ 86). 

The difference between the calculated one-neutron 
binding energies of the fragment at scission B* and 
ground state B n is plotted in Fig. [131 We find that this 
difference can be as large as 2 MeV in absolute value. 
Such differences are a consequence of the evolution of 
single-particle neutron gaps as a function of deformation. 



2. Neutron multiplicity 

The multiplicities calculated from Eq. ([8j arc shown in 
Figs. [14] - US as functions of fragment mass for the four 
studied nuclei. On figures [T4l and \W\ a solid line has been 
added to guide the eye. Typical saw-tooth structures are 
observed, displaying maxima and minima. These struc- 
tures appear correlated with the quadrupole deformation 
of fragments at scission. 

In the case of 226 Th fission, the neutron multiplicity 
curve displays pronounced structures separated by five 
mass units from Af iag = 110 to Af rag = 150, that are 
linked to: i) fragment deformations (Figs. 151 and ITD|) . ii) 
the deformability of the fragments (the softness of po- 
tential energies with respect to axial quadrupole and oc- 
tupolc deformations), and iii) the one- neutron binding 
energy (see Fig. [T2")) . 

For Fm isotopes, the curves look more regular, and 
show that neutron emission is almost vanishing around 
^4-frag = 130 and is maximum around Af lag = 120. In 
Fig. 1151 comparison is made with the experimental data 
for spontaneous fission of 256 Fm (4?| ■ The agreement be- 
tween theoretical values and measurements is rather sat- 
isfactory, as the global data pattern is well reproduced. 
However, calculations appear to underestimate the num- 
ber of emitted neutrons in the ^4f rag = 90 — 130 region 
and a second minimum is found around j4f rag = 144. As a 
consequence, the calculated number of emitted neutrons 
is 30% smaller than experimental data. 

These discrepancies probably come from our model as- 
sumptions. Part of the overall underestimation of the 
number of emitted neutrons is presumably due to the 
fact that the deformation energy has been calculated with 
constraints placed only on axial quadrupole and octupolc 
deformations of fragments. More realistic calculations 
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FIG. 14: (Color online) 226 Th. Calculated neutron multiplic- 
ity as a function of fragment mass. The solid line is to guide 
the eye. 
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FIG. 15: 256 Fm. Neutron multiplicity versus fragment mass. 
Comparison between predictions (solid symbols) and data [47t ] 
(empty symbols). 



should include the effect of higher order multipole frag- 
ment deformations such as q^ and q 60 . As for the second 
minimum at A{ rag ~ 144, we think it may originate from 
the method adopted to define the two fragments, espe- 
cially for cases where many particles are present in the 
neck. Handling the neck rupture in a more realistic man- 
ner and including higher order multipole deformations in 
the calculation of deformation energies would change our 
predictions. However we have doubt that these consider- 
ations, alone, would be pertinent enough for bringing in 
significant improvements. These discussion will be fur- 
ther extended below. 
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FIG. 16: (Color online) Neutron multiplicities for 258 Fm (up- 
per panel) and 260 Fm (lower panel). Lines are to guide the 
eye. 



G. Deviation from the unchanged charge 
distribution 



Introduced in 1962 by Wahl, the fragment unchanged 
charge distribution Ztjcd (i-e- charge polarization) is the 
charge number of a fragment with a given mass .Af ra g, if 
its Z/A ratio were the same as the one of the fissioning 
nucleus (Hi: 



^UCD 



Z{ s .Af Ia 
AZ 



(9) 



The deviation AZf rag = Zf rag — Ztjcd of the charge 
of the fragment Zf rSug from the unchanged charge distri- 
bution Ztjcd is plotted in Figs. [T7] and [TH] as a function 
of the fragment charge for 226 Th and 258 Fm. Values of 
AZf rag for 256 Fm and 260 Fm are very close to the ones of 
258 Fm. We first observe that AZf lag is globally positive 
for light fragments and negative for heavy ones. This fca- 
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FIG. 17: Nascent fission fragment proton pairing energy (up- 
per panel) and deviation from unchanged charge distribution 
(bottom panel) as functions of fragment charge for 226 Th. 



ture stems from the fact that heavy systems may sustain 
stronger neutron excess than light ones, as observed and 
discussed for several fissioning systems ^| Hi} ■ The pat- 
terns displayed by AZf lag as functions of Zf rag are quite 
contrasted. While both sets of AZf lag values show sharp 
structures as Zf rag increases, AZf rag globally decreases 
in 226 Th from AZ frag ~ 1 to AZ frag ~ -1. In con- 
trast, in 258 Fm the AZf rag values reach a plateau with 
AZf rag | ~ 0.5 for Zf rag > 57 and Zf Tas < 43. In an at- 
tempt to understand the origins of these sharp structures 
and different global trends, we have sought for possible 
correlations with other structure properties, namely pro- 
ton separation energies of: i) fissioning nuclei along scis- 
sion lines, and ii) each nascent fragment. No clear cut 
correlation is found. However the structures observed in 
the AZf rag values for both nuclides seem to coincide with 
the variations of the fragment pairing energies Ep an . ( p ) 
as can be seen when comparing the plots in the upper 
and bottom panels in Figs. [TTl and [T51 



H. Total kinetic energy 

1. Distance between fragments 

The total kinetic energy (TKE) of a given fragmenta- 
tion can be estimated from the formula 



TKE 



e 2 Z n Z^ 
d c h 



(10) 



where e is the electron charge, Zh (Z\) the charge of 
the heavy (light) fragment, and d c h the distance between 
fragment centers of charge at scission. The distance d c h 
deduced from our calculations is plotted as a function of 
fragment mass for 226 Th, 256 Fm, 258 Fm and 260 Fm in 
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FIG. 18: Same as Fig. [17] for 258 Fm. 
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FIG. 19: Distances between nascent fragment centers of 
charge calculated as functions of fragment mass for 226 Th and 

256-260p m ^ 



Fig. [TSJ For all considered nuclei, the distance between 
fragment centers of charge at scission falls in the range 
dch = 14 — 20 fm. The distance d cm between centers of 
mass has also been calculated. The difference Sd between 
d cnl an d dch appear rather small: 8d ~ 0.08 fm in 226 Th 
and Sd ~ 0.05 fm in 256-260 Fm _ 



2. Total kinetic energy for 226 Th 

The TKE values of 226 Th fission fragments are plot- 
ted as functions of fragment mass in Fig. [201 The dots 
represent the result of Eq. (fT0)l whereas the solid line 
follows the experimental data of Ref. [4(| HH obtained 
in electro-magnetic induced fission measurements. One 
notices that theoretical results present many more struc- 
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FIG. 20: 226 Th. TKE values of nascent fission fragments as a 
function of fragment mass. Comparison between predictions 
(solid dots) and data [13, EH]. 



tures than do experimental data. This difference may be 
explained from the fact that the experimental measure- 
ments correspond to an excitation energy of the fissioning 
nucleus of the order of 11 MeV, whereas formula (flC)]) is 
valid only for low energy fission. As well known, an in- 
crease in the fission energy smooths out kinetic energy 
distribution. In particular the kinetic energy in the sym- 
metric mass region increases [52j which explains why ex- 
perimental TKE display only a very shallow minimum 
for A frag = A/2. 

The experimental TKE values display maxima for 
A bag ~ 132 (Z frag ~ 54) and A frag ~ 94 (Z frag ~ 36), 
whose positions coincide with those found in our calcula- 
tions. The sharp structures in the theoretical TKE stem 
from the compactness of scission configurations marked 
with the symbols d, e and f in Fig. [5] They do not 
show up in the experimental TKE values, as details of 
the energy landscape along the scission line are presum- 
ably washed out by the dynamics of the Coulomb induced 
fission process (53|. Nonetheless, theoretical results are in 
qualitative agreement with experimental data, with de- 
viations never exceeding 15%, and the calculated mean 
TKE value (TKE) th ~ 1 69 Me V, is found close to the ex- 
perimental mean value (TKE) cxp = 167.7 ±3.4 MeV 0]. 



3. Total kinetic energy for 



'Fro. 



The TKE values of 256 - 26 0p m g ss i on fragments are 
plotted in Fig. [21] as functions of fragment mass. The 
TKE curves look rather similar in all three isotopes. 
They display a sharp peak reaching TKE ~ 250 MeV 
for symmetric fission. These features are characteristic 
of compact scission, where the fissioning system gives rise 
to nearly spherical fragments separated by a small dis- 
tance. Furthermore, we also observe that the full width 



FIG. 21: (Color online) 256 " 260 Fm. TKE values of nascent 
fission fragments as functions of fragment mass. 
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FIG. 22: 256 Fm. TKE values of fission fragments as func- 
tions of fragment mass. Results of theoretical calculations 
(dots) are displayed with pre-neutron-emission data (contour 
diagram). The solid line represents the experimental average 
TKE Hj. 



Atke at half maximum is narrowing from Atke = 20 u 
to Atke = 14 u in going from 256 Fm to 260 Fm. 

Theoretical TKE in 256 Fm (black points) arc compared 
with experimental ones [54| in Fig. [55] Calculated results 
are found in very good agreement with the average TKE 
data for asymmetric fission (Af rag = 138 — 150). However, 
they overestimate the experimental values for Af rag < 138 
by up to 16%. This latter feature may indicate that the 
calculated distance between fragment charges c? c h is too 
small for scission configurations close to symmetry. This 
remark is consistent with the fact that for the same frag- 
mentations, ^f rag values calculated from deformation en- 
ergies is underestimated (sec Fig. 1151) . These under- and 
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over-estimations are interpreted as due to our study, re- 
stricted to the (920 1 930 ) deformations, which favors com- 
pact scission. Accessing elongated fission configurations 
for nearly-symmetric fragmentations of 256 Fm implies 
that at least three collective coordinates must be con- 
sidered in constrained HFB calculations. 



IV. CONCLUSION 

In this work, large scale HFB calculations using the 
Gogny D1S force have been performed in order to in- 
vestigate structure properties of 226 Th and 256_260 Fm at 
scission and the characteristics of fission fragments along 
scission lines. Scission configurations are first analyzed 
assuming that axial quadrupole and octupole collective 
coordinates play a major role in fission. We have found 
from our constrained HFB calculations that the scission 
mechanism depends on which heavy nuclide and frag- 
mentation are considered. This mechanism may display 
either a smooth or an abrupt character in the (9201930) 
plane. The former property means that the potential en- 
ergy of the fissioning system changes smoothly over de- 
formation from outer saddle to scission and beyond where 
Coulomb repulsion takes place between fission fragments. 
This scission property is found for the Fm symmetric 
fragmentations. For asymmetric fragmentations in all 
nuclei of present interest, there is a sudden drop in po- 
tential energies whenever scission takes place. To accom- 
modate with these contrasted properties, post-scission 
points arc defined for matter densities present in the neck 
that are weaker than p = 0.06 fm~ 3 . With this criterion, 
scission lines are calculated, and the fragmentations de- 
termined assuming sharp cuts across the necks. 

Properties of fission fragments and correlations with 
properties of fissioning systems along scission lines have 
been discussed. These comprise potential and deforma- 



tion energies, quadrupole and octupole deformations, to- 
tal kinetic energies, prompt neutron emissions, deviation 
from unchanged charge distribution, and energy parti- 
tioning. All these properties reflect either shell and/or 
pairing contents of potential energies of both fission frag- 
ments and fissioning nuclei, in particular for multipolc de- 
formations, neutron multiplicities, and total kinetic ener- 
gies. Predictions are found in reasonably good agreement 
with experimental data for total kinetic energy ( 226 Th, 
256 Fm) and prompt neutron multiplicity ( 256 Fm) of fis- 
sion fragments. 

The present microscopic analysis shows that the struc- 
ture of the two-dimensional (q2o, 930) potential energy 
surface in 226 Th is similar to those previously calculated 
in U and Pu. The different behavior with respect to scis- 
sion found in Fm isotopes and the fact that symmetric 
elongated configurations do not appear in our description 
may indicate that a collective space with more that two 
dimensions is needed to describe scission configurations 
and fragment properties in these nuclei. Preliminary in- 
vestigations show that other heavy actinides probably 
also require an enlargement of the dimension of the col- 
lective space used. In view of the encouraging results 
obtained so far, in particular in light actinides, it seems 
worth attempting to extend the present static calcula- 
tion to three dimensions or even more. Of course, a 
description of fission observables such as fragment mass 
and kinetic energy distributions will also require to ex- 
tend the two-dimensional dynamical model employed in 
Uranium [2(j to higher dimensions. 
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